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We describe how the finite-difference time-domain (FDTD) technique can be used to compute the quasinormal 
mode (QNM) for metallic nano-resonators, which is important for describing and understanding light-matter 
interactions in nanoplasmonics. We use the QNM to model the enhanced spontaneous emission rate for dipole 
emitters near a gold nanorod dimer structure using a newly developed QNM expansion technique. Significant 
enhanced photon emission factors of around 1500 are obtained with large output /3-factors of about 60%. (c) 
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Resonant cavity structures have the ability to trap light 
in very small spatial volumes which has a wide range 
of applications in nanophotonics [1]. Various miniatur¬ 
ized cavity structures have been developed over the years 
to manipulate light at the subwavelength scale, and ex¬ 
treme nanoscale confinement is now possible with metal¬ 
lic nano-resonators (MNRs). In a frequency regime near 
a localized surface plasmon (LSP) [2], the local density 
of optical states (LDOS) can be increased dramatically. 
Consequently, the spontaneous emission (SE) rate of a 
dipole emitter can be enhanced via the Purcell effect. 
The resonant enhancement from MNRs has applications 
in chemical sensing [3], high resolution imaging [4,5], 
optical antennas [6] and single photon emission [7]. 

While the optical properties of MNRs are being ac¬ 
tively pursued, numerically modelling of the basic cav¬ 
ity physics is extremely demanding, and analytical solu¬ 
tions of the modes only exists for very simple structures 
such as spheres. For resonant cavity structures, the nat¬ 
ural modes of the system are called quasinormal modes 
(QNMs) [8,9], defined as the frequency domain solutions 
to the wave equation with open boundary conditions (the 
Silver-Miiller radiation condition). Kristensen et al [10] 
first used the QNMs to introduce a rigorous definition 
of the “generalized effective mode volume” and Purcell 
factor [11], and applied these results to photonic cavity 
structures. For MNRs, the QNMs also form the natural 
starting point for developing analytical theories of light- 
matter interactions in nanoplasmonics [12-14]. 

One of the most common numerical techniques for 
obtaining the cavity mode for dielectric cavities is the 
finite-difference time-domain (FDTD) technique [15]. 
The FDTD technique allows one to simulate open 
boundary conditions with “perfectly matched layers” 
(PMLs) located at the leaky mode region outside the 
cavity. For dielectric cavities, this open-boundary FDTD 
approach has been shown to yield excellent agreement 
with direct integral equation methods [10]. Other time- 
domain techniques such as the discontinuous Galerkin 


time-domain approach can also use PMLs [16]. For met¬ 
als, a major problem occurs when using the usual mode 
calculation approach with a dipole excitation source, 
since the extracted mode depends sensitively on the 
dipole position [17], and is therefore incorrect. Thus it 
has been common practice to excite the MNR with a 
plane wave source and obtain the scattered field. How¬ 
ever, this scattered field is not the same field as the QNM 
and it cannot be properly normalized for use in quantum 
optics, e.g., for obtaining the Purcell factor and effective 
mode volume [12]—two well known quantities that help 
describe the underlying physics of cavity light-matter 
interactions. While some frequency-domain techniques 
exist for computing the QNMs of MNRs [18,19], it is 
highly desirable to be able to compute the QNMs using 
the commonly employed and general FDTD technique. 

The FDTD method is already widely used by the plas- 
monics community, and its accuracy for obtaining the 
enhanced field has been verified against other numeri¬ 
cal techniques such as the multipole expansion technique 
[20]. In addition, the LDOS can be calculated by employ¬ 
ing a dipole excitation source [21-23], which can also 
model local field effects, e.g., associated with finite-size 
photon emitters inside a MNR [21]. While direct dipole 
calculations are feasible, they are very time consuming 
and require a new simulation (which may require many 
hours of computational time) for each spatial position of 
the dipole emitter; thus a QNM picture would be much 
more efficient, since it allows one to simulate dipole re¬ 
sponses both as a function of space and frequency as 
soon as the QNM is obtained and properly normalized. 

In this Letter we first describe how the FDTD tech¬ 
nique [24] can be efficiently employed to obtain QNMs of 
a MNR by filtering the scattered field with a temporal 
window function. We compute the spatial dependence 
of the QNM and the effective mode volume, and show 
excellent agreement with full-dipole calculations for the 
enhanced SE factor of a dipole emitter. We then show 
how a gold nanorod dimer can act as an efficient single 
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photon source with large Purcell factors (1500) and im¬ 
pressive output /3-factors (around 60%). In contrast to 
spherical dimer structures [25], we find that the nanorod 
dimer acts to increase the /3-factor for good photon emit¬ 
ters (in comparison to a single resonator), and has the 
additional advantage of yielding resonant frequencies in 
the visible spectrum. 

The QNM has a complex eigenfrequency, Cb^ = — 

iy^, with a quality factor Q = The QNM is 

normalized through [8,9] 
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Since the eigenfrequencies are complex, the QNMs di¬ 
verge in space and each part of Eq. (1) diverges, but 
the total sum converges quickly in space [12,14]. This 
convergence occurs approximately when the outgoing 
field becomes purely oscillatory, rather than evanescent. 
For spatial points near the resonator, the transverse 
part of the photon Green function can be expanded 
as [9,14] GT(ri,r 2 ;w) = One 

can then derive the effective mode volume and en¬ 
hanced SE factor Fa(ra,co’), where is the spatial po¬ 
sition of a dipole emitter. For dipole positions near the 
resonator, these quantities are defined through [10,14] 


and 
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Purcell factor, Fp, multiplied by a factor r] to account 
for any deviations at from the field maximum ro, cav¬ 
ity polarization, and cavity resonant frequency [14]; here 
Ua is the unit vector in the direction of the dipole emit¬ 
ter, ub is the background-medium refractive index in 
which the MNR is embedded and Ac is the corresponding 
wavelength. Equation (3) assumes a single QNM expan¬ 
sion, which is valid for the MNR light-matter interaction 
regimes that we study below. 

The FDTD technique can also obtain the enhanced 
SE factor directly, but only at the dipole position, which 
can be used to check the numerical accuracy of the QNM 
equations. Specifically, one can define [26] 


pFDTDr , r _ Im{na-G™™(r„,ra;c^)Tia} 

lm{na • GB(ra,ra;cj) • 11^} 

where Gb is the known homogeneous medium Green 
function [26], and G™™(ra, r^; cj) is obtained from 
FDTD at the dipole source position [22]. Note that our 
full-dipole FDTD calculations have been shown else¬ 
where to yield good agreement with analytical results 
for metal spheres [27] (within 5% agreement over broad¬ 


band frequencies and various spatial positions). 

Next we discuss how FDTD is usually used to obtain 
resonator cavity modes and highlight the main problem 
for metallic structures. As a direct space-time simulation 
method, FDTD can obtain the spectral response of sys¬ 
tem from a Fourier transform of its time response. Thus 
for any incident field with a finite bandwidth, the com¬ 
puted scattered field, E^(cc;), will contain the resonant re¬ 
sponses of the system, e.g., spectral peaks corresponding 
to the QMNs. Defining the total electric field as E^°^(co’) 
(including the incident source), then the scattered field 
is simply E^(co’) = E^°^(co’) - E^(cl;), where E^(co’) is the 
homogeneous solution without the MNR (i.e., no scatte¬ 
ring). For a typical dielectric cavity with a high Q res¬ 
onance, almost any incident field or dipole source field 
can be used to obtain the mode of interested if one ap¬ 
plies a temporal window to subtract off the excitation 
source. A common method of choice is to use a dipole 
excitation source and begin the Fourier transform of the 
scattered field after the excitation pulse has gone, which 
has been shown to yield good agreement with other mode 
solving techniques [10]. However, for a metal, the dipole 
source does not efficiently excite the QNM because of 
large losses, and extra care is needed in applying a tem¬ 
poral window function to compute the QNM of interest. 
Below we describe how the FDTD technique can obtain 
the true QNM, if the QNM is spectrally well isolated. 

To cover a broad spectral range, we use a plane-wave 
excitation source with a 6 fs time duration. To efficiently 
excite the QNM, we choose the polarization of the source 
field to be parallel to the expected polarization of the 
QNM (e.g., along the rod axis for a nanorod). After ob¬ 
taining the scattered field, we compute the QNM from 
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where oOc is the real part of the eigenfrequency 
chc, 2Vln2 Twin gives the FWHM (full-width at half¬ 
maximum) of the time window, and toff is the time offset 
from the center of the source pulse to the center of the 
time windowing function. One criterion for the selection 
of time window parameters is to set cjc^win/^^off ^ 
while keeping Twin as large as possible. To understand 
why this works, consider the ideal case for which the inte¬ 
gral in Eq. (5) is carried out from —oo to oo, and assume 
the scattered field is given by E^(t) = 

((<^i,7i) % if ^ % 4)5 after performing the integra¬ 

tion, the amplitude of the Xth component is proportional 

^ (i^ ^ ^off "^2 

to jE^e 2 |e with = uji — and 

— 7i— ^off/^win* resouauces of the system are 

well separated, so that \Auji\r^iri ^ 1, and lAcj^l ^ lAy^j 
for i % c, then the QNM resonance, Cbc = — iTc will 

be the only surviving term after applying the windowing 
function. This is achieved by having toff/r^in ~ 

_ *off 

the QNM of interest; however, due to the factor e ^win , 
toff should not be too large, otherwise the spectrum will 
be too weak and influenced by numerical errors. We have 
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Fig. 1. Spatial profile of the dipole mode of a gold 
nanorod dimer, showing both the scattered field and 
QNM. Left: |E^(x, 0; cjc)!, and right: |f (x, 0; cjc)! at 

ujc = 291.06 THz. The excitation source is a ^-polarized 
plane wave incident in the x direction. 

found that a good choice is to simply choose toff = 47r/7c, 
i.e., two times the lifespan of the QNM. 

Next we apply the above technique to investigate a 
gold dimer structure made up of two identical nanorods. 
We choose a rod radius r^ = 15 nm with an axis 
length I = 100 nm (along y). We use the Drude model, 

s{uj) = 1 — with parameters similar to gold, 

with cjp = 1.26 X 10^^ rad/s (plasmon frequency), and 
7 = 1.41 X 10^^ rad/s (collision rate). In order to get a 
larger enhancement of the SE for a quantum dipole, we 
consider two nanorods that are parallel with each other. 
The separation gap is set to 20 nm, which helps mini¬ 
mizes nonradiative decay and also allows sufficient space 
to embed a quantum emitter such as a quantum dot (see 
dipole arrow in Fig. 2). From FDTD analysis, we find 
the dipole mode is around a;c/27r = 291.06 — i20.28 THz 
(cjc ~ 1.2 eV), which is redshifted (by about 34 THz) 
with respect to a single nanorod [14] due to the bonding 
effect; the corresponding quality factor is Q ~ 7.2, which 
is smaller than that for a single nanorod (Qsingie ~ 10). 

To obtain the scattered field and the QNM, a y- 
polarized plane wave with frequency ci;/27r = 291.06 THz 
is employed, incident in the x direction. The simu¬ 
lation domain size is 12 x 8 x 6 (Fig. 1) and 

2.4 X 2.4 X 2.4 ym^ (Fig. 2), and we use a conformal 
meshing scheme with a maximum step size 40 nm in all 
directions; a smaller refined mesh of 1 nm is used around 
the metal dimer; 100 layers of PML have been used with 
symmetric (antisymmetric) boundary condition in z {y) 
direction for both simulation and the time step is 1.8875 
as. We use a time window with parameters toff = 100 fs 
and Twin = 23.3/Vln2 fs. The left and right panels of 
Fig. 1 show the scattered field and QNM, respectively. 
As anticipated, the QNM is seen to have an increasing 
field value for positions further away from the resonator, 
which is caused by the outgoing boundary conditions and 
the complex eigenfrequency. A close-up view of the QNM 
near the resonator is shown in the left panel of Fig. 2, 
and the correspondent mode volume is calculated to be 
Feff ~ 1.9 X 10~^{Xc/nB)^ (which is about double that of 
the single nanorod). The oscillating charge distribution 
is of opposite sign at both ends of each nanorod; also, 
the oscillating charge at the bottom of the upper nanorod 
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Fig. 2. Close up view of the QNM profile of the gold 
dimer. Left: |f (x, 0; cjc)!, middle: |fy(x, 0; cjc)!, and 

right: \fx{x,y,0]UJc)\ at uJc = 291.06 THz. A ^-polarized 
dipole at the center of the gap is shown by the white dot 
and double arrow in the left panel. 


has a different sign from the charge at the top of the lower 
nanorod, which makes the localized surface plasmons of 
the individual nanorods couple effectively. Consequently, 
the field inside the cavity (gap) is significantly enhanced 
and ^-polarized. The middle and right panels in Fig. 2 
show the y and x components of the QNM at 2 : = 0. The 
node of the ^-component sits around both ends of the 
nanorods, slightly inside this metal. Similar nodal lines 
appear for an electric dipole composed of charge =Fg. 

We remark that for these calculations it is important 
to choose a ^-polarized plane wave, which efficiently ex¬ 
cites the QNM of interest. If we use an x-polarized plane 
wave, with the parameters same as above, then a rather 
strange pseudo-mode is obtained. The failure of using 
x-polarized plane wave to obtain the QNM can be ex¬ 
plained by the fact that the QNM can not be efficiently 
excited since the dipole moment lies along the y direc¬ 
tion. So some care is needed in choosing the correct exci¬ 
tation source, though in practice this is easy to do when 
exciting the QNM of interest. We have also checked that 
our mode calculation approach successfully reproduces 
the mode profile and QNMs for dielectric cavity struc¬ 
tures, e.g., in Ref. 10 and also for other MNRs. 

Single photon enhancement factors. As a possible ap¬ 
plication of using MNRs for single photon emitters, a 
large enhancement of the SE factor is desired. This can 
be achieved by having a quantum dipole at the gap center 
of the dimer. In Fig. 3, the enhanced SE factor, of a y- 
polarized dipole emitter at the center of the gap is shown 
by the blue (dashed) curve using a full-dipole numerical 
calculation with no approximations, i.e., F™™ [22,23]. 
Next we use the QNM to obtain the enhanced SE fac¬ 
tor [14]. The result is shown by the orange (solid) curve 
in Fig. 3; it agrees extremely well with the full-dipole 
FDTD calculations. Indeed, the agreement also confirms 
that the entire response is dominated by a single QNM 
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Fig. 3. Enhanced SE factor for the gold nanorod dimer 
structure, with the blue (dashed) curve given by direct 
EDTD calculation (Eq. (4)) and the orange (solid) curve 
given by the mode expansion technique (Eq. (3)). The 
emission dipole is ^-polarized at the center of the gap of 
the dimer (see left panel of Eig. 2.) 


which is highly desired for single photon source appli¬ 
cations. Erom our results, we obtain impressive peak 
emission factors of around 1500. Note that although the 
QjV factor of the dimer structure is about three times 
smaller than the single nanorod, the emission factor of 
the dimer is around three times larger than for a single 
nanorod at an equivalent dipole position. Thus the Q/V 
is clearly not the main metric to describe the enhanced 
emission factors for dipole away from the field antin¬ 
ode; note the reason for not choosing the field antinode, 
is that the emission would be severely quenched from 
Ohmic heating and a single mode Purcell factor ceases 
to have any meaning [14]. In this regard, the other im¬ 
portant figure-of-merit for emitting single photons is the 
output-coupling p factor (probability of photon emission 
via radiative decay in the far field) of the dipole emit¬ 
ter, which is computed here to be /3 ~ 58%; this again 
is better than the single gold nanorod, by about 10%; 
with respect to single nanorod, the electric field of the 
dimer is repelled outside the lossy nanorods (especially 
around the hot spots) due to opposite charge distribution 
between them. This leads an overall nonradiative decay 
reduction (oc J Im[e(r, cj)] |E^°^(r, co’)pdE), even though 
the volume of the integral is doubled. 

In summary, we have described how the commonly 
employed EDTD method can be used to efficiently obtain 
the QNMs for metallic resonators by applying a filtering 
function to the scattered field. These calculations are 
verified by comparing the resulting SE emission factor 
with full-dipole EDTD numerical calculations, where we 
find excellent agreement. Using this technique, we have 
proposed a metal nanorod dimer structure that acts an 
efficient single photon source, yielding a large Purcell 
factor of 1500 and an output [3 factor of around 60%. 
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